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Quark model results for the B —* tt, p decays are analysed, making use of the dispersion formulation 
of the model: The form factors at q 2 > are expressed as relativistic invariant double spectral represen- 
tation over invariant masses of the initial and final mesons through their light-cone wave functions. The 
dependence of the results on the quark model parameters is studied. For various versions of the quark 
model the ranges T(B° -> vr+ez^) = (7±2) x 10 12 |K&| 2 a -1 , T(B° -» p + eu)/T(B° -> vr+eP) = 1.45±0.1, 

\& ■ and r^/r^ = 0.7 ± 0.08 are found. The effects of the constituent quark transition form factor are 

^ , briefly discussed. 

2 ! PACS number(s): 13.20.-v, 12.39.Hg, 12.39.Ki. 

i*r ■ Weak decays of hadrons provide an important source of information on the parameters of the 

^ standard model of electroweak interactions, the structure of weak currents, and internal structure of 

' hadrons. Hadron decay rates involve both the Cabibbo-Kobayashi-Maskawa matrix elements and 
hadron form factors, therefore the extraction of the standard model parameters from the experiments 

\ on hadron decays requires reliable information on hadron structure. 

' A theoretical study of hadronic matrix elements of the weak currents inevitably encounters the 

problem of describing the hadron structure and requires a nonperturbative consideration. This gives 
. the main uncertainty to the theoretical predictions for hadron transition amplitudes. 

For the transitions between hadrons each containing a heavy quark and light degreess of freedom, 
the number of the independent form factors is considerably reduced due to the heavy quark symmetry 
[Q]. For instance, in the leading 1/mg order the transition between heavy mesons is described in terms 
of the single Isgur-Wise function, which should be estimated within a nonperturbative approach. The 
O • 0(1/ nig) corrections to this picture can be consistently calculated within the Heavy Quark Effective 
Theory Q, an effective theory based on QCD in the limit of large quark masses (a detailed review 
D ! can be found in ||]). 

For the decays caused by the heavy-to-light quark transitions the argumentation of the heavy-to- 
heavy transitions does not work, and the situation turns out to be much less definite. For instance, 

* i 

for the decays B — > tt, p the uncertainty of the results of nonperturbative approaches such as the 



quark model QCD sum rules [ 10 1— [ 22] , and lattice calculations pql is too large to draw 



any definite conclusion on the form factor values and decay rates (see Tables |], |2|) 

A step forward in the understanding of heavy-to-light transitions was recently done by B.Stech 
who noticed that relations between heavy-to-light form factors can be obtained if use is made of 



the constituent quark picture [16]. Namely, assuming that (i) the momentum distribution of the 
constituent quarks inside a meson is strongly peaked with a width corresponding to the confine- 
ment scale, and (ii) the process in which the spectator retains its spin and momentum dominates 
the transition amplitude, [] the 6 form factors can be reduced to a single function just as it is in the 
case of the heavy-to-heavy transition. These relations are expected to be valid up to the correc- 
tions O (2m u MB / (Mg + M 2 — q 2 ))- Although these corrections cannot be estimated numerically, the 
relations can be a guideline to the analysis of the heavy-to-light decay processes. 

Constituent quark picture has been extensively applied to the description of the decay processes 



[p.7| . Although in the first models by Wirbel, Stech, and Bauer (WSB) || and Isgur, Scora, 



1 Actually, one more assumption on the dynamics of the procees is employed. Namely, the picture in includes 
both the constituent and current quarks. And for deriving the final relations it is important that the momentum of 
the current quark coincides with the momentum of the corresponding constituent quark. This assumption allows one to 
avoid the appearence of the constituent quark transition form factor which should be taken into account if the picture 
with only constituent quarks is considered. 
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Grinstein and Wise (ISGW) || quark spins were not treated relativistically, it has become clear soon 
that for a consistent application of quark models to electroweak decays, a relativistic treatment of 
quark spins is necessary ||. The exact solution to this complicated dynamical problem is not known, 
but a simplified self-consistent relativistic treatment of the quark spins can be performed within the 



light-cone formalism [18] . The only difficulty with this approach is that the applicability of the model is 
restricted by the condition q 2 < 0, while the physical region for hadron decays is < q 2 < (Mi — Mf) 2 , 
M{ j being the initial and final hadron mass, respectively. So, for obtaining the form factors in the 
physical region and decay rates and lepton distributions, assumptions on the form factor g 2 -behavior 



were necessary. A procedure to remedy this difficulty has been proposed in [19]. 



The approach of []T9| is based on the dispersion formulation of the light -cone quark model [20|. 
Namely, the transition form factors obtained within the light -cone formalism at q 2 < [17], are 
represented as dispersion integrals over initial and final hadron masses. The transition form factors 
at q 2 > are derived by performing the analytic continuation in q 2 from the region q 2 < 0. As a 
result, for a decay caused by the weak transition of the quark Q(rrii) — > Q(rrif), form factors in the 
region q 2 < (m,j — nif) 2 are expressed through the light-cone wave functions of the initial and final 
hadrons. We apply this approach to the analysis of the B — > (ir, p) decays and study the dependence 
of the results on the quark model parameters, such as constituent quark masses and wave functions. 
We also check the fulfillment of the Stech relations in particular model calculations. 

The amplitudes of the semileptonic decays of a pseudoscalar meson P(M\) into the final pseu- 
doscalar P(M 2 ) and vector V(M 2 ) mesons have the following structure |l]] 

< P(M 2 ,p 2 )\V fl (0)\P(M l ,p 1 ) > = f + (q 2 )(pi+p 2 )^ + f-(q 2 )(pi-p 2 )^ 
<V(M 2 ,p 2 ,e)\V ll (0)\P(M l , Pl )> = 2g(q 2 )e^e* v {p 2 )p^ 2 (1) 
<V(M 2 , P2 ,e)\A^0)\P(M 1 , Pl ) > = if(q 2 )e^ + ia + (q 2 )(6*p 1 )(p 1 +p 2 ) l , + ia^q 2 )(€*p 1 )(p 1 - P2 )^ 

We denote both the pseudoscalar and vector meson masses as M 2 but the relevant value is taken in 
each case. 

Our goal is to calculate the form factors of the transition between S-wave mesons at < q 2 < 
(Mi — M 2 ) 2 within the constituent quark picture. In this picture the initial meson P(M\) is a bound 
state of the constituent quarks Q{m 2 )Q{m 3 ), the final meson is an S-wave bound state Q{mi)Q{m 3 ), 
and the transition process is described by the graph of Fig.l. 

We however start with the region q 2 < and make use of the fact that the transition form factors 
calculated within the light-cone quark model [|l7|] can be written as double spectral representations 
[|l9|].[^0| over the invariant masses of the initial and final mesons 

oo s+(s 2 ,q 2 ) 

f ,2\_t , 2\ f ds 2 G 2 {s 2 ) f cisiGi(si) fi{si,s 2 ,q z ) 
/l(<?) - /2l(9) J tt( S2 - Mf) J vr( Sl -M 2 ) 16AV2( Sl>a2l(? 2)' (2) 

(mi+m 3 ) 2 a - (s 2 ,g 2 ) 



where 



sf(s 2 ,q 2 



2ml 



w f ( 2 2 2\ 2/ 2 i 2\ , I 2 2\/ 2 2\ i \l/2/ 2 2\\l/2/ 2 2 2\\ 

x ( S 2{Q — mi — m 2 ) — q {m l + m 3 ) + [m l — m 2 )[m l — m 3 ) ± A ; {s 2 , m 3 , m 1 )\ 1 [q ,m 1 ,m 2 )J 
and 

A(si, s 2 , s 3 ) = (si + s 2 - s 3 ) 2 - 4sis 2 - 

Here G± 2 are the vertex functions which describe the constituent quark structure of the mesons, 
/2i(? 2 ) is the form factor of the constituent quark weak transition m 2 — ► mi. In what follows we 
adopt the conventional approximation f 2 \{q 2 ) = 1. The same representations of the form factors has 



been obtained in [21] taking into account the contribution of two-particle singularities in the Feynman 
graphs. 
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The double spectral densities fa of the form factors f + ,g,a+, and /, which give a nonvanishing 
contribution to the cross-section of the semileptonic decay in the case of zero lepton mass, read 
(hereafter S3 = q 2 ) 



f + (s 1 ,s 2 ,s 3 ) = D + (ai + a 2 )D 3 



9{Sl,S 2 ,S 3 ) 



m\OL 2 + m 2 ai + 01.3(1 — ol\ — a 2 ) 



2ft 



^ + m 1 + m 3 

a + (si, s 2 , s 3 ) = 2 [2m 2 an + 2m 3 (ai - a n ) - mia 2 - m 2 (ai - 2ai 2 ) 

Cat + 63(011 + a<i2 



/(si,s 2 ,s 3 ) 



-m 3 (l - ai - a 2 - 2a i2 ) + 
M 2 



/d(si, s 2 , s 3 ) + M 2 a + (si, s 2 , s 3 ) 



s 2 + mi + m 3 

si-s 2 -s 3 M 2 -M^-s 3 



/d(si,s 2 ,s 3 ) = 4 



£ + 2/3(m 2 - m 3 ) + 



2^ 

C 3 /3 



2M 2 



^/si + mi + m 3 



where 



L> = si + s 2 - (m 2 - m 3 ) 2 - (mi - m 3 ) 2 , D 3 = s 3 - (mi - m 2 ) 2 - L>, 
C = si + s 2 - (m 2 - m 3 ) 2 - (mi + m 3 ) 2 , C 3 = s 3 - (mi + m 2 ) 2 - C, 

m 2 . 2 2\ ^1 / 2 2\ ^3 / 2 2\ 

E = mim 2 m 3 + -^-(s 2 - m x - m 3 ) + -y («! ~ ™, 2 - m 3) ~ ~y ( S3 ~ m i ~ m 2)> 



1 



«i 



a 2 



Hsi,s 2 ,s 3 ) 
1 

A(si,s 2 ,s 3 ) 
1 r 



(si + s 2 - S3)(s2 - m 2 + m|) - 2s2(si - m 2 , + m|" 
(si + s 2 - s 3 )(si - + m|) - 2si(s 2 - m 2 + m|" 



/3 = - 2m| — ai(si — m| + mj) — a 2 (s 2 — m 2 + m|) 



an = a 2 + 4/3- 



•52 



ai 2 = aia 2 — 2ft 



si + s 2 - s 3 



(3) 
(4) 
(5) 

(6) 



A(si,s 2 ,s 3 )' " A(si,s 2 ,s 3 )' 

For an S-wave meson (vector and pseudoscalar) with the mass M built up of the constituent quarks 
mi and m 2 , the function G is normalized as follows [19| 



G 2 (s)ds A 1 / 2 (s,m 2 ,m; 



2j ( S -(mi-m 2 ) 2 ) = 1. 



(7) 

tt(s-M 2 ) 2 8tts v v 7 ' K ' 

Notice that the double dispersion representation without subtractions are valid for the form factors 
/+, /_ g, a + , and a_, while the form factor / requires subtractions. 

Now, we obtain the form factors at q 2 > by performing the analytic continuation in q 2 . For the 
function which has at q 2 = s 3 < the structure 



<Ks 3 ) 



ds 2 G 2 (s 2 ) f dsiGi(si) 



(m 1 +mz)' 1 



vr(s 2 - Mf ) J vr(si - M{ 



Po(si,s 2 ,s 3 ) | Pijsi, s 2 ,s 3 ) | P 2 (si,s 2 ,s 3 ) 



AV2(si,s 2 ,s 3 ) A3/2( Sl , S2) s 3 ) AV2( S1 , S2)S3 ) 



(8) 

where Pi are polynomials of s, the analytical continuation to the region q 2 > yields the following 
expression at q 2 < (m 2 — mi) 2 



«3 



ds 2 G 2 (s 2 ) f dsiGi(si 



(mi+m-i) 2 



vr(s 2 - Mf) J vr(si - Mf ) 



Pq(si,s 2 ,s 3 ) | Pi{si,s 2 ,s 3 ) | P 2 (si,s 2 ,s 3 ) 



A 1/2 (si,s 2 ,s 3 ) A 3 / 2 (si,s 2 ,s 3 ) A 5 / 2 (si,s 2 ,s 3 ) 



0) 



3 



-26(s 3 ) 



ds 2 G 2 (s 2 ) _ 
vr(s 2 - Mf) J tt 



si 



My/2 



I / \ , ^i(si)-fo(sf) , 2 (si) - 2 Of ) - 0' 2 («f - sf )" 
&o(,siJ H 5 H 



si 



,-R\2 



s 3 + mf 



2 



+ 



\ 



S3 + m\ 



rrio 



+ (mf - mf), s 3 < (m 2 - mi)' 



(y^ + \/si) 2 > A(si, s 2 , s 3 ) = (si - sf )(si - s 
Gi(s)P„(si,s 2 ,s 3 )6 , (si < sj~) 



f), and 



4>n(s) 



(Si -M2)(si -sf)™ +1 /2 



(10) 



(11) 



For deriving this expression it is important that the functions (?i 2 have no singularities in the r.h.s. 
of the complex s— plane [21]. Along with the normal Landau-type contribution connected with the 
subprocess when all intermediate particles go on mass shell (the first term in @), the anomalous 
contribution (the second term in (||)) emerges at q 2 > |22]. The normal contribution dominates the 
form factor at small positive q 2 and vanishes at the 'quark zero recoil point' q 2 = (m 2 — mi) 2 . The 
anomalous contribution is negligible at small q 2 and grows as q 2 — > (m 2 — mi) 2 . 

We are now in a position to apply these results to the B — > ir,p decays, in which case m 2 = m^, 
mi = 7713 = m u . To this end the quark model parameters (m u ,mb and the vertex functions) should 
be specified. 

For a pseudoscalar meson built up of quarks with the masses mi and m 2 , it is convenient to 
introduce the function w related to the vertex function G as follows 



GOO 



7T 



mt 



m. 



2\2 



M 1 



\f2 \Js- (mi - m 2 )' z 



,3/4 



-w(k), 



A 1 / 2 (s, mf , m 2 ,) 



The normalization condition (^) for G yields the following normalization condition for 



w 2 (k)k 2 dk = 1. 



(/.' 



(12) 



(13) 



The function w is the ground-state 5-wave radial wave function of a pseudoscalar and vector meson 
for which a simple exponential form is usually chosen 



w(k) = exp(-fc 2 /2/3 2 



(14) 



We consider several sets of the quark model parameters used for the description of the meson spectra 
and elastic form factors (Table ||) . 

For analysing the results of calculations, introduce the functions Ri(q 2 ) such that 



U(Q 2 ) 
V(q 2 ) 

Mq 2 ) 

A 2 (q 2 ) 



R+(q 2 ), 

(Mi + M 2 )g(q 2 ) 



1 



[l+r)R v (q 2 ), 
l + r 2 



y 



(Mi+M 2 ) J v * ' l + r 
-(Mi+M 2 )a + (q 2 ) = (l + r) 



1-r 2 



(l+r) 



^R2(q 2 ) 



where r = M 2 /Mi, y = q 2 /M 2 . As found by Stech fHf , Ri(q 2 ) should be equal up to the corrections of 
the order O (2m^Mi/(M 2 + Mf — q 2 )). For the decay B — > it the corrections about 10% are expected 
even at q 2 = 0. 

Table |] presents the parameters of the best fits to the calculated form factors in the form 

Ri(0) 



R; 



(i - q 2 /m 2 y 
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These fits approximate the form factors with better than 0.5% accuracy in the relevant kinematic 
range and can be used for the calculation of the decay rates. In all the form factors, the anomalous 
contribution comes into the game only at q 2 > 15 Gev 2 ; below this point it is negligible. The 
decay rates are calculated with the form factors from Table || via the formulas from Ref. J?]] . Table |5| 
summarises the results on the form factors and decay rates. 
The following conclusions can be drawn: 

(i) Relativistic effects are important in the B —* tt, p decays, as our results considerably differ from 
those of the models with nonrelativistic treating quark spins with the same parameters (Sets 1 and 
2). 

(ii) The functions Ri are equal with an expected 10-20% accuracy as found by Stech for any set of 
the parameters. However, the magnitude and g 2 -dependence of the functions R strongly depend on 
these parameters (see Fig. 2). The 'quark model average values' for the decay rates can be determined 

as 

T(B° -► ir+ev) = (7 ± 2) x W 12 \V ub \ 2 a" 1 , 
F(B° -» p+eP)/F(B° -> -K+ev) = 1.45 ± 0.1, 

r L /r T = 0.7 ±0.08. (15) 

These values are not far from the predictions of other models (Table ^), except for the results by 
Narison ( |i~2j| ). This is mainly due to a very specific decreasing q 2 -behavior of A\ in ||i~2|| . 

The ratio of the decay rates have rather good accuracy, while their absolute values are less definite. 
Perhaps, the way out of this situation lies in a more detailed consideration of the constituent quark 
transition form factor. We use an approximation /21 = 1. However, there are arguments in favour of a 
nontrivial g 2 -dependence of /21. Actually, within the constituent quark picture, the constituent quark 
form factor arises from the two sources. First, it appears as a bare form factor which describes the 
constituent quark amplitude of the weak current defined through current quarks. One can assume this 
bare form factor to be close to unity at relevant q 2 . Secondly, this bare form factor is renormalised 
by the constituent quark rescatterings (final state interactions) ^0|. These very interactions yield 
meson formation and are not small at least in the region of q 2 near the meson mass. So the relation 
/21 = const cannot be valid for all q 2 . Whether it is strongly violated or not at the q 2 of interest is not 
clear yet. Setting /21 = 1 in @ means that we take into account the effects of the constituent quark 
binding in the initial and final meson channels, and neglect the quark rescatterings in the g 2 -channel. 
The constituent quark form factor should depend on the parameters of the model. One can expect 
that once the constituent quark form factor is taken into account thoroughly, the results on the meson 
transition form factors and decay rates will be weakly dependent on the quark model parameters. In 
the ratio of the decay rates the constituent quark form factors cancel, and these results seem to be 
more reliable. 

I am grateful to V. Anisovich, S. Godfrey, I. Grach, I. Narodetskii, and K. Ter-Martirosyan for 
valuable discussions and to V. Nikonov for help in numerical calculations. The work was supported 
by the Russian Foundation for Basic Research under grant 95-02-04808a. 



2 The decay rates are rather sensitive to the details of the g 2 -behaviour of Ai. In the light-cone approach the form 
factor Ai has a specific feature: unlike other form factors it cannot be determined uniquely. Using the dispersion language, 
the form factor Ai requires subtractions which cannot be fixed uniquely. This is a general situation for any amplitude 
with a vector particle in the initial or final state. Such an amplitude is transverse with respect to the momentum of the 
vector particle, and this yields a necessity of subtractions in some of the form factors. However, different definitions of 
the subtraction procedure lead to rather close results. 
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Table 1: The form factors of the decays B — > ir, p at q 2 = 0. The labels QM, SR, and LAT stand for 



Quark Model, Sum Rules, and Lattice, respectively. 



Ref. 


/+(0) 


V(0) 


Ai(0) 


A 2 (0) 


^(o)Mi(o) 


A 2 (0)/Ai(0) 


QM 


WSB m 
ISGW § 


0.33 
0.09 


0.33 
0.27 


0.28 
0.05 


0.28 
0.02 


1.2 
5.4 


1.0 

0.4 


SR 


BBD |1C] 
Ball (111 
Narisongl 


0.24± 0.025 
0.26± 0.02 
0.23± 0.02 


0.6±0.2 
0.45±0.05 


0.5±0.1 
0.38±0.04 


0.4±0.2 
0.45±0.05 


l.lli 0.01 


l.lli 0.01 


Lat 


APE |lja 
APE [3b 
ELC Ha 
ELC Qb 


0.29± 0.06 
0.35± 0.08 
0.26± 0.16 
0.30± 0.19 


0.45±0.22 
0.53±0.31 
0.34±0.10 
0.37±0.11 


0.29±0.16 
0.24±0.12 
0.25±0.06 
0.22±0.05 


0.24± 0.56 
0.27± 0.80 
0.38± 0.22 
0.49± 0.26 


2.0±0.9 
2.6±1.9 
1.4±0.2 
1.6±0.3 


0.8±1.5 
1.0±3.1 
1.5±0.7 
2.3±0.9 



Table 2: Decay rates of the transitions B — > (tt + , p + )e v in units |Vub| 2 x 10 12 s 



Ref. 


T(B -> vr) 


T(B p) 




Tl/Tt 


WSB [| 


7.43 


26.1 


3.5 


1.34 


ISGW § 


2.1 


8.3 


3.95 


0.75 


Ball PI 


5.1± 1.1 


12±4 




0.06±0.02 


Narison|l|] 


3.6± 0.6 




0.9±0.2 


0.2±0.01 


APE [13| 


8±4 









Table 3: Parameters of the quark model. 



Parameter Set 


m u 


m b 




Pb 


Setl i 


0.35 


4.9 


0.4 


0.4 


Set 2 § 


0.33 


5.12 


0.31 


0.41 


Set 3 fL7j 


0.29 


4.9 


0.29 


0.386 


Set 4 [23] 


0.25 


4.7 


0.32 


0.55 



Ta 


ble 4: 


Parameters of the best fits to the ca 


culated form factors. 


Set 


R+(0) 


M+ 


n + 


R v (0) 


My 


ny 


fli(0) 


Mi 


m 


R 2 (0) 


M 2 


n2 


Set 1 


0.29 


6.29 


2.35 


0.30 


6.28 


2.36 


0.27 


7.07 


2.65 


0.25 


6.13 


2.17 


Set 2 


0.20 


6.22 


2.45 


0.20 


6.22 


2.46 


0.20 


6.78 


2.65 


0.19 


6.00 


2.34 


Set 3 


0.21 


5.90 


2.33 


0.21 


5.90 


2.35 


0.21 


6.50 


2.70 


0.20 


5.90 


2.45 


Set 4 


0.26 


5.44 


1.72 


0.29 


5.46 


1.73 


0.29 


5.68 


1.67 


0.28 


5.36 


1.67 



Table 5: Form factors of the B° — > (tt + , p + )e v decays at q 2 = and the calculated decay rates in 
units \V ub \ 2 x 10 12 8~K 



Set 


/+(0) 


V(0) 


Mo) 


A 2 (0) 


T(B -» vr) 


T(B - p) 






Set 1 


0.29 


0.34 


0.24 


0.21 


9.2 


12.6 


1.38 


0.77 


Set 2 


0.20 


0.22 


0.17 


0.16 


5.0 


7.25 


1.45 


0.73 


Set 3 


0.20 


0.24 


0.19 


0.17 


6.1 


8.8 


1.44 


0.62 


Set 4 


0.26 


0.33 


0.26 


0.24 


8.9 


13.8 


1.55 


0.67 
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Figure 1: One- loop graph for a meson decay. 



O 
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10 20 n 2 



Figure 2: The g 2 -dependence of the form factor / + for the decay b — > ir at various quark model 
parameters. Set 1 - upper dashed line, set 2 - lower solid line, set 3 - dash-dotted line, set 4 - lower 
dashed line, the monopole formula with M po i e = Mb* = 5.34 GeV - dotted line. 
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